Neutron wave packet tomography 
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A tomographic technique is introduced in order to determine the quantum state of the center of 
mass motion of neutrons. An experiment is proposed and numerically analyzed. 
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In experimental physics one often faces the following 
question: "Given the outcomes of a particular set of mea- 
surements, which quantum state do they imply?" Such 
inverse problems may arise for instance when setting up 
and calibrating laboratory sources of quantum states, or 
in the analysis of decoherence and other deteriorating 
effects of the environment, or in some special tasks in 
quantum information processing such as eavesdropping 
on a quantum channel in quantum cryptography. 

The dctermination of the quantum state represents a 
highly nontrivial problem, whosc history can be traced 
back to the early days of quantum mechanics, namely to 
the Pauli problem Q, Q ; the experimental validation had 
to wait until quantum òptics opened a new era. The the- 
oretical predictions of Vogel and Risken 0] were closely 
followed by the experimental realization of the suggested 
algorithm by Smithey et al. Since then many im- 

provements and new techniques have been proposed: an 
up-to-date overview can be found in Ref. 0. Recent 
progress in instrumentation has made it possible to ap- 
ply these techniques to a variety of different quantum 
systems such as fields in optical cavities, polarization and 
external degrees of freedom of photons, or motional states 
of àtoms. 

In this Letter we propose an experiment for deter- 
mining the quantum states of the center-of-mass motion 
of neutrons. In accordance with quantum theory, these 
massive partides can be associated with a wave function 
describing their motional state. Neutrons are suitable 
objeets for many quantum mechanieal experiments due 
to their interaction with all four bàsic forces, the ease 
of detecting them with almost 100% emeieney, and their 
small coupling to the environment 0|. In marked con- 
trast with light, neutron vacuum field and thermal back- 
ground can usually be ignored. This makes it possible, 
for instance, to prepare superpositions of macroscopically 
separated quantum states — the so-called Schròdinger cat 
states — that would be extremely dimcult to realize with 
other quantum systems because of their fragility with re- 
spect to decoherence. In all experiments performed so 
far, the existence of the Schròdinger cat states of neu- 
trons has been indircctly demonstrated via interferomet- 
ric effects, but the full evidence for the nonclassicality of 
these states, including the presence of the negative vàlues 



of the reconstructed Wigner function, is still missing. 

In the following, we will first briefly review the present 
neutron interferometric techniques and the means of cre- 
ating highly nonclassical motional states of neutrons. In 
the second part of the Letter, an experiment will be pro- 
posed for the complete reconstruction of these quantum 
states. 

Neutron tomography - The set of measurements that 
can be done on neutrons to determine their quantum 
state is severely limited by the very low time resolution of 
the available detectors. In quantum òptics, this obstacle 
can be overcome by mixing the weak input field with a 
strong local oscillator. By changing the phase <f> of the 
oscillator one can measure the spectral decompositions 
of all quadratures, 
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x and p being the canonically conjugated operators of 
position and momentum. Of course, no such local oscil- 
lators exist for neutrons. However, notice that massive 
partides experience a transformation of the type Q in 
the course of free evolution: x(t) = x + (p/m)t, where 
to is the mass. Thus free evolution of the wave packet 
followed by a position sensitive measurement yields infor- 
mation about a subset of quadratures X^,, <f> 6 [0, 7r/2]. 
Free evolution was utilized e.g. for the reconstruction of 
transversal motional states of helium àtoms in ■ Here 
we are interested in the longitudinal degrees of freedom. 
Since neutron detectors have very bad time resolution, 
free evolution alone cannot be used to generate a tomo- 
graphically complete set of measurements. 

Feasible measurements on thermal neutrons consist of 
measurements of the contrast and phase of interference 
fringes in an interferometric setup, see Fig. ^ (without 
momentum kick) , and also spectral analysis of the neu- 
tron beam using an adjustable Bragg-reflecting crystal 
plate together with a position sensitive detector. This 
set of observables is not tomographically complete be- 
cause the measurable (complex) contrast of the interfer- 
ence pattern @ (H = 1), 



T(Ax) = (^|e íAa:í5 |^) = / \a(p)\ 2 e l/lxp dp, (2) 
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is not sensitive to the phase of a(p) = (p\ip), and no 
information about quadratures other than p is available. 
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FIG. 1: Scheme of a perfect crystal neutron interferometer. 
The incoming beam is split at the first crystal plate, reflected 
at the middle plate and recombined again at the third plate. 
The detector is placed in beam O where the visibility is higher 
due to the same number of reflections/transmissions. In ad- 
dition to a position shift Aa; routinely used in neutron exper- 
iments, a momentum kick Ap has been added in path II in 
order to make the interferometric measurement tomographi- 
cally complete; see text. 



Obviously, the situation would be different if one could 
shift both the position (phase) and the momentum of the 
incoming wave packet inside the interferometer. Such 
a thought experiment is shown in Fig. ^ I n that case, 
the Wigner function describing the ensemble of measured 
neutrons would be related to the measured contrast 



r(Ap, Ax) = tr {pé Apx é Axp } 

e l{Ap)x {x\p\x + Ax)dx 

by a simple integral transformation, 

W(x,p) -- 



(3) 



e- l — +mx+ivp T(-u,v)dudv, (4) 



where u = Ap, v = Ax and p is the state to be recon- 
structed. 

Although this thought experiment looks simple, its 
experimental realization, according to Fig. ^ would be 
rather difficult. Large momentum kicks acquired by the 
neutron in the lower arm would change its de Broglic 
wavelength and spoil the Bragg reflection at the last 
crystal plate. Therefore we will now propose a modi- 
fied scheme that can substitute the interferometric setup 
of Fig. H| 

Setup - In the new scheme shown in Fig. the incom- 
ing neutrons polarized in the +z direction, |W) = |t/>)|z + ), 
where ip denotes the spatial degrees of freedom, first 
propagate freely through a distance L, undergoing a uni- 
tary operation Ui = exp[— ip 2 L/(2p )]. In the following, 
we will assume that the input wave packets are quasi- 
monochromatic, <r p /po ~ a few percent, with spread a p 
and central momentum pq. This condition guarantees 
that the action of the RF coil is practically eqivalent to 
a momentum "kick." The generalization to more general 
states will be presented elsewhere. After the region of 
free propagation, the neutrons are let through an RF 
coil placed in a static magnètic field polarized along the 
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FIG. 2: Setup for the tomography of motional states of neu- 
trons. Mi, M2 - magnètic mirrors; L - region of free propaga- 
tion; B - static magnètic field controlling the momentum kick 
Ap; F - box containing a static magnètic field (aligned along 
+y) and an RF coil; D - detector; the state is "prepared" 
by P, which can be a chopper or a monochromator. The ar- 
rows denote the polarization of neutrons after reflections from 
magnètic mirrors. 



— y direction, see Fig. [3 As a result of the interaction 
between the neutron and the coil, the y+ compo- 
nent of the input state will be decelerated (accelerated) . 
Assuming that the region of interaction is short, so that 
the dispersion of the wave packet of the neutron can be 
neglected in the coil, in the quasi- monochromatic approx- 
imation, the net momentum transfer can be described by 
the effective unitary operator, 

U-2 = e- iA ^ 2 \y.)(y + \ + e iA ^' 2 \y + ){y.\, (5) 



where 
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Prior to detection, the partides are polarized along the 
+z direction again, so as to erase the which-way infor- 
mation stored in the polarization degree of freedom. The 
probability of a neutron being detected is given by the 
norm of the transmitted component, 

P = Tr{n(Ap,Ax)p}, (7) 

where p refers only to the spatial degrees of freedom and 

Tl(A P ,Ax) = (z + \uM\z + )(z + \U 2 U 1 \z + ) (8) 

= (l + e l ^ p2 e lApx e~ l ^ p2 )/4 + h.c. 

= (l + e lAp(x+L·p/po) )/'í + h.c. (9) 

= (1 + e iAxí> e iAp£ e- iAxAp / 2 )/4: + h.c., (10) 



where we denoted 



Aa; = 



ApL 2p,BmL 
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Dircct inversion - The detection probability (Q reads 
P(Ap, Ax) = i + ÍRe{r(Ap, Ax)e 2lAxA P}. (12) 
Since the beam is quasi-monochromatic one has for Sx — 

T(Ap, Ax + Sx) ~ r(Ap, Ai)e™/ 2 (13) 

from which the imaginary part of the complex degree of 
coherence T(Ap, Ax) can be obtained. 

Summarizing, the tomography of a neutron state con- 
sists in the following four steps: 

(i) A set of pairs of independent variables {Bj,Lj} is 
chosen covering a certain range B € [0, B max ] and L 6 

[0, imax] ■ 

(ii) For each pair Bj , Lj the shifts Apj in (JüJ) and Axj 
in 111(1 are calculated, and the corresponding intensities 
P(Apj, Axj) are measured with and without an auxiliary 
shift Sx = n/2. 

(iii) The complex degree of coherence T(Apj, Axj) is cal- 
culated from the two intensities using Eqs. (|12|) - (|13fl . 

(iv) Finally, the Wigner function of the input neutrons is 
calculated with the help of inversion formula where 
the integrals are approximated by sums over Axj and 
A Pj . 

According to l@J, the contrast T(Ap, Ax) is essentially 
the Fourier transform of the Wigner function W(x,p). 
Therefore the largest vàlues of Ap and A:r are related 
to the smallest resolved details in x and p respectively. 
Namely (reinserting h), 

ApMAX = h/Sx min , AlMAX = Tl/Spmm, (14) 

where Sx min and Sp mm denote the x and p resolutions. 
By © and ijTTJl one gets 
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Spnún = "F^n 



(15) 



For a neutron of wavelength Ao = 0.37nm assuming 
the reasonable vàlues Lmax = lm and -Bmax = 0.1T one 
gets Sx m i n — 60/im and Sp m i n = U x 10 m . 

Radon inversion - It is interesting to give an alternative 
interpretation of the proposed measurement in Fig. [21 
Notice, that the POVM elements in Eq. JÜJ can also be 
restated in terms of quadrature operators, 

U(A Pl Ax) = (1/4) (l + e*"* 8 ) +h.c, (16) 
where (in hxed units) 

Ax L 

Xg = cos 9x + sin 6p, tan 6* = — — = — , (17) 

Ap po 



the quadrature Xg, 

P(Ap,Ax) = l/2 + Re{C X9 (w)}/2, (18) 
{CxM) = I Px e (x)e^dx. (19) 



and u> — y/ Ax 2 + Ap 2 . Thus, for a fixed 9, the data 
contain information about the characteristic function of 



By changing L one changes the quadrature measured, 
whilc lo, which depends on both L and B, determines the 
observed spatial frequeney of the probability distribution 
of this quadrature. The observed quadratures range from 
x (for L = 0) to p (for L — > oo). From the measurement 
of Cx e ( w ) ; the "shadows" Px g (x) of the Wigner function 
can be obtained by the Fourier transformation, which 
in turn yield the Wigner function by an inverse Radon 
transformation. This is an alternative way of reconstruc- 
tion the Wigner function from the measured data in the 
setup Fig. [3 

Statistical inversion - The procedures outlined above, 
based on the direct inversion formula have several 
drawbacks: (i) Realistic data are always noisy. In that 
case, formula (@J can yield unphysical results, such as 
the Wigner representation of a non-positive defmite op- 
erator. (ii) The Wigner function in Eq. J3J depends on 
the measured data indirectly, through the complex degree 
of coherence T, which itself has to be estimated with the 
help of an auxiliary position shifter. This intermediate 
step is, certainly, not necessary as all available informa- 
tion about the Wigner function of the incoming neutrons 
is contained in the raw data measured without any aux- 
iliary position shift. In order to avoid these problems, we 
propose to use the maximum-likclihood quantum state 
reconstruction [5J, |9j, llOfl . The main advantages of this 
method compared to the above direct inversion are: (i) 
Asymptotically, for large data samples it provides the 
best performance available. (ii) Any prior information 
about the measured neutrons and the known statistics of 
the experiment can be used to increase the aceuracy of 
the reconstruction. (iii) The existing physical constraints 
can be easily incorporated into the reconstruction. Most 
notably, this technique guarantees the positivity of the 
reconstructed density operator. (iv) It can be applied 
directly to raw counted data. 

Assuming that the statistics of the experiment is Pois- 
sonian, the maximum-likclihood reconstruction amounts 
to minimizing the Kullback-Lcibler distance (relative en- 
tropy) between the measured data /(Ax, Ap) and the 
renormalized theoretical probabilities p(Ax, Ap) / J^p of 
Eq. (|12J) ■ As has been shown in |3, llfj , the maximum- 
likely density matrix can be obtained as a hxed point of 
the iterations of a nonlinear operator map. 

As follows from the parameter estimates given after 
Eq. (TSJ, the proposed tomography scheme using ther- 
mal neutrons will likely have sufhcient resolution in mo- 
mentum. On the other hand, even for well monochro- 
matized thermal beams, the resolution in position is ex- 
pected to be worse (possibly even by several orders of 
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FIG. 3: Reconstructed Wigner functions of Gaussian states 
(upper row) and superpositions of Gaussian states ("cats;" 
lower row) from the simulated data. A 50 x 50 matrix of 
Ax and Ap shifts was used for the maximum-likelihood in- 
version. Left column: reconstructed original states; middle 
column: reconstructed time-evolved states with a resolution 
5s m i n /i co h = 1/2; right column: reconstructed time-evolved 
states with reduced resolution, <5x m in/ícoh = 10. Grey is the 
zero level, white (black) represents positive (negative) vàlues. 



magnitude) than typical coherence lengths. The simu- 
lations in Figure [31 illustrate the effect of the restricted 
range of Ap on the reconstruction. Consider first the re- 
construction of a minimum uncertainty Gaussian wave 
packet in its moving frame, parameterized by its co- 
herence length Z coll , I^g) ck / exp(-fc 2 ^ oh )|fc)cífc. (The 
choice of a minimum uncertainty state is only for illus- 
trative purposes.) Provided the apparatus has a suffi- 
cient spatial resolution, Sx m i n < l co h, a faithful recon- 
struction is readily obtained, see the upper left panel. 
More realistic measurement with Sx m i n > Z co h would ob- 
viously yield a Wigner function smoothed out along the x 
axis. However, the states measured in a real experiment 
are not going to be minimum uncertainty states. The 
experimenter will rather deal with time evolved states 

*(T)} oc J exp{ik 2 T/2m-k 2 l 2 oh )\k)dk that are strongly 
affected by dispersion. As a consequence, the wave packet 
spread very soon becomes larger than the resolution 
limit, Sxt ~ T/(ml C oh) ^> <5x m i n , and a good reconstruc- 
tion can be achieved with a realistic apparatus. Compare 
the upper middle and right panels, showing reconstruc- 
tions with a sufficient resolution fo m j n = Z C oh/2 and a 
reduced (but more realistic) resolution Sx m i n = 10/ C oh- 

The imaging of non-classical states is a much more 
delicate problem. Let us consider the superpositions of 
spatially separated Gaussian states (Schròdinger cats), 

*cat) oc [1 + exp(ipA)]|í , G). Such states can be pre- 
pared e.g. by means of a double loop perfect crystal in- 
terferometer . As has been shown |é| , the preparation 
of thermal neutron cat states is possible, with separations 
exceeding the corresponding coherence lengths of the in- 



dividual components A 3> ^coii- Provided the apparatus 
has sufficient resolution, the nonclassicality of this state 
is manifested by the negative regions of the reconstructed 
Wigner function, see the lower left panel of Fig. Tak- 
ing dispersion into account, the ordering of the relevant 
parameters is Z co h < A < 5xt, and one can easily obtain 
femin < Sxt- As the simulations show, a realistic mea- 
surement whose position resolution is much worse than 
the coherence lengths of the individual cat state compo- 
nents tends to wipe out the negative regions of the re- 
constructed Wigner function; compare the lower middle 
and right panels of Fig. [3] On the other hand, the main 
features of such exòtic states, such as their non-Gaussian 
character and also the global spatial properties of which 
little is known today, should still be accessible to a re- 
alistic wave packet tomography. To resolve more subtle 
quantum interference effects of the order of the coher- 
ence length, more refined experimental techniques may 
however be needed. An idea could be to replace thermal 
neutrons by ultracold neutrons, for which much larger 
momentum shifts Ap (and thus much smaller Sx m - m , pos- 
sibly even smaller than A) can be obtained. 

In conclusion, wc have proposed and analyzed an ex- 
perimental scheme for determining the motional states 
of neutrons. With the help of a magnètic field and free 
propagation, this apparatus realizes quadrature measure- 
ments on neutrons by measuring overlaps of the two 
transformed components of the initial state. This is an 
analog of the quantum optical homodyne detection in 
neutron òptics, achieved without the use of a strong co- 
herent source of neutrons. 
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